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ABSTRACT 

We present uniformly measured star formation histories (SFHs) of 60 nearby {D < 4 Mpc) dwarf 
galaxies based on color-magnitude diagrams of resolved stellar populations from images taken with 
Hubble Space Telescope and analyzed as part of the ACS Nearby Galaxy Survey Treasury program 
(ANGST). This volume-Hmited sample contains 12 dSph/dE, 5 dwarf spiral, 28 dirr, 12 dSph/dlrr 
(transition), and 3 tidal dwarf galaxies. The sample spans a range of ^ 10 in Mb and covers a wide 
range of environments, from highly interacting to truly isolated. From the best fit SFHs we find 
three significant results: (1) the average dwarf galaxy formed > 50% of its stars by z ^ 2 and 60% 
of its stars by z 1, regardless of current morphological type; (2) the mean SFHs of dis, dTrans, 
and dSphs are similar over most of cosmic time, and only begin to diverge a few Gyr ago, with the 
clearest differences between the three appearing during the most recent 1 Gyr; and (3) the SFHs are 
complex and the mean values are inconsistent with simple SFH models, e.g., single bursts, constant 
SFRs, or smooth, exponentially declining SFRs. The mean SFHs are in general agreement with the 
cosmic SFH, although we observe offsets at intermediate times (z ~ 1) that could be evidence that 
low mass systems experienced delayed star formation relative to more massive galaxies. The sample 
shows a strong density-morphology relationship, such that the dSphs in the sample are less isolated 
than dIs. We find that the transition from a gas- rich to gas-poor galaxy cannot be solely due to 
internal mechanisms such as stellar feedback, and instead is likely the result of external mechanisms, 
e.g., ram pressure and tidal stripping and tidal forces. The average transition dwarf galaxy is slightly 
less isolated and less gas-rich than the typical dwarf irregular. Further, the transition dwarfs can be 
divided into two groups: interacting and gas-poor or isolated and gas-rich, suggesting two possible 
evolutionary pathways. 

Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: formation — galaxies: stellar 
content — color-magnitude diagrams (HR diagram) 



1. INTRODUCTION 

Dwarf galaxies have come to play an increasingly im- 
portant role in understanding how galaxies form and 
evolve. As the smallest, least luminous, and most com- 
mon systems in the universe, dwarf galaxies span a wide 
range in physical characte ristics and occupy a diverse 
set of environment s (e.g.. iZwic kv' 1957'; Hodge 1971'; 
Koo fc Kroiil 119921: [Marzke fc da Costa 1997: Mateq 
19981: Ivan den Berghl 120001 : iKarachentsev et all I2004D . 
making them excellent laboratories for direct studies of 
cause and effect in galaxy evolution. The low average 
masses and metallicities of dwarf galaxies suggest they 
may be the best available analogs to the seeds of hierar- 
chical galaxy formation in the early universe. 

A cohesive picture of the evolution of dwarf galax- 
ies remains elusive. Historically, evolutionary scenar- 
ios have often been considered in the context of a dual 
morphological classification, namely dwarf spheroidals 
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(dSphs; we include dwarf ellipticals in this general cat- 
egory) and dwarf irregulars (dIs). The former are clas- 
sified based on a smooth morphology, with no observed 
knots of SF, and are generally found to be gas-poor. The 
latter exhibit morphological evidence for current /recent 
SF activity and have a high gas fraction. A third, 
rarer class of so-called 'transition' dwarf galaxies, have 
a high g as fraction yet little to no recent SF activ- 
ity (e.g., Sandage & Hoffman 1991; 'Skillman & Bended 
1995: Mateo 1998; M iller et al . 2001; Dolphin et al. 2003 
Young et al.( l2007riDellenbusch et all |2008). These 
galaxies ma y be an evolutiona ry link between dIs and 
dSphs (e.g.. iGrebel et al.1l2003| ) or simply could be ordi- 
nary dIs witnessed between massive star forming events 
(e.g., Skillman et al. 2003a). Understanding the rela- 
tionship between these three types of dwarf galaxies, 
specifically determining if and how dIs evolve into dSphs, 
is among the n iost pressing questions in dwarf galaxy 
evolution (e.g., 'Baade' 1951^ 'Hodgel fT97ll; 'Dekel & SiM 
[1986; Binggcfi 1986; Hodge 1989; Skillman & Bcndep 
fl995; Mateo 1998; Dckcl & Woo 2003; Grcbcl et al. 2003 
i Rico t ti fc Gnedin 2005; May er et al. 2006; Orban et at] 
l2008i: iKazantzidis eTaLll2010t ). 

Resolved stellar populations have proven to be an 
incredibly powerful tool for observationally constraining 
scenarios of dwarf galaxy evolution. Past patterns of 
SF and chemical evolution are encoded in a galaxy's 
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optical color-magnitude diagram (CMD). To extract 
this informat ion, a number of sophisticated algo- 
rithms (e.g., 'Tosi et al.' '19891: iTolst oy fc Sahal 119961: 

lel T Il997t iH ohzman et aP 
2001; 



Gallart ct al. 1996; Mighel 
1999t iHernandez e t al 



DolDhid I200l 



20071: lYukfcLig 



1999; 
Ikuta & Arirnoto. 
2007t 



1997r 



Harris & Zaritskv 

, , 20021: ICo lc ctaL 

. lAparicio fc Hidalgo! l200"a 

Cignoni &: Tosill20l0 ) have been developed to measure 



the star formation history (SFH), i.e., the star forma- 
tion rate (SFR) as a function of time and metallicity, 
by comparing observed CMDs with those generated 
from models of stellar evolution. The robustness of 
this method has become increasingly solidified with 
a variety of techniques and stellar models converg- 
ing on consistent solutions for a ran ge of galaxies 
(e-g., .S killman & Gallart 2002; SkiUm an et all l2003bl : 
iGallart et al.l l2005: Monclh ct al.i i2010a| b^ 

Analysis of CMD-based SFHs have become particu- 
larly prevalent in studies of the formation and evolution 
of the Local Group (LG). The LG contains ~ 80 dwarf 
galaxies ranging froni highly isolated t o str ongly interact- 
ing (e .g.. lMateolll998tl^an den Berga2000l : lTolstov et al.l 
Hooa). Results from CMD-based SFHs of dwarf galaxies 
in the LG have revealed that both dSphs and dis feature 
complex SFHs, typically with dominant stell ar compo- 
nents older than 10 Gyr (e.g.,|Ho dgc 1989; Matedll998l: 
IDolphin et aT][200l iTolstov et al, i2009.1 . The SFHs of 
individual LG dwarf galaxies, as well as aggregate com- 
pilations, now serve as the basis for our understand- 
ing of the evolution of dwarf galaxies and have signifi- 
cantly ad vanced our knowled ge of galactic group dynam- 
ics fe.g- IMayer et al.' '200Tal.lbl, 120061; 'Or ban et all 120081: 
iTolstov et a'il2b 09: Mavcr 2010; Kazantzi dis et al.lboiOD . 

Beyond the LG, only a small subset of dwarf galax- 
ies have expl i citly measured CMD-ba sed SFHs (e.g.. 
Dolphin c FaD 120031: iWeisz et al.l 120081 : iMcQuinn etahl 
2009i , ,2010i) . Ground-based observations of resolved stel- 



lar populations in more distant galaxies are challenging, 
due to the faintness of individual stars and the effects 
of photometric crowding. However, observations from 
the Hubble Space Telescope (HS T), particula rly the Ad- 
vanced Camera for Surveys (ACS: lFord erall[l998l' l . have 
revolutionized the field of resolved stellar populations, 
producing stunning CMDs of dwarf galaxies both in and 
beyond the LG. Outside the LG, past HST observations 
of resolved stellar populations have been fairly piecemeal, 
with individual or small sets of galaxies as the typical 
targets. Subsequent analysis often employed different 
methodologies or sought different science goals. This lack 
of uniformity also makes it challenging to place the re- 
sults from LG studies in the context of the broader uni- 
verse, as LG dwarf galaxies may not be r epresentative of 
the la rger dwarf galaxy population (e.g.. Ivan den Berghl 
[2003). 

The ACS Nearby G alaxy Survey Treasury (ANGST; 
iDalcanton et al.ll2009[ ) was designed to help remedy this 
situation. Using a combination of new and archival imag- 
ing taken with the HS T/ACS and Wide Field Plane- 
tary Camera 2 (WFPC2: lH"oltzman et al. 1995), ANGST 
provides a uniformly reduced, multi-color photometric 
database of the resolved stellar populations of a volume- 
limited sample of nearby galaxies (D < 4 Mpc) that are 
strictly outside the LG. Of the ^ 70 galaxies in this sam- 
ple, 60 are dwarf galaxies that span a range of ~ 10 in 



Mb and that reside in both isolated field and strongly 
interacting group settings, providing an unbiased statis- 
tical sample in which to study the detailed properties 
dwarf galaxy formation and evolution. 

In this paper, we present the uniformly analyzed SFHs 
of 60 dwarf galaxies based on observations, photometry, 
and artificial star tests produced by the ANGST pro- 
gram. The focus of this study is on the lifetime SFHs of 
the sample galaxies, with the recent (< 1 Gyr) SFHs th e 
subject of a separate paper (|D. Weisz et al..llin Prep.l ). 
We first briefly review the sample selection, observations, 
and photometry in Sj2l We then summarize the technique 
of measuring SFHs in fJH] In 521 we discuss and compare 
the resultant SFHs in the context of both dwarf galaxy 
formation and evolution. We then explore our results 
with respect to the morphology-density relationship in 
Sj5l Finally, we discuss the evolution of dwarf galaxies in 
the context of cosmology in [JH Cosmological parameters 
used in this paper assume a sta i idard WMAP-7 cosmol- 
ogy as detailed in iJarosik et al.l (j2010[ ). 

2. THE DATA 

In this section, we briefly summarize the selection of 
the ANGST dwarf galaxy sample along with the obser- 
vations and photometry. A more detailed discus sion of 
the A NGST program can be found in iDalcanton et al.l 
(l2009l) . 

2.1. Selection and Final Sample 

We constructed the initial list of ANGST galax- 
ies based on the Catalo g of Neighboring Galaxies 
(jKarachentsev et al.l I2004D . consisting exclusively of 
galaxies located beyond the zero velocity surface of the 
LG ([van den Bcrgh 2000). We selected potential tar- 
gets with |6| > 20° to avoid observational difficulties 
associated with low Galactic latitudes. By simulating 
CMDs and crowding limits, we found that a maximum 
distance of ^ 3.5 Mpc provides the optimal balance be- 
tween observational efficiency and achieving the program 
science goals. Because the sample of galaxies within 
3.5 Mpc contains predominantly field galaxies, we chose 
to extend the distance limits in the direction of the 
M81 Group. This extension added to the diversity of 
galaxies in the sample, while maintaining the goal of 
observational efficiency, due to the M81 Group's close 
proximity to the fiducial distance limit (Dmsi ~ 3.6 
Mpc) and low foreground extinction values. Similarly, 
we also included galaxies in the direction of the NGC 
253 clump (DAf253 ~ 3.9 M pc) in the Sculptor Filament 
(|Karachentsev et al.ll2003[) . further extending the range 
of environments probed by the sample, while maintaining 
strict volume limits for galaxy selection. 

In this paper, we analyze a sample of 60 ANGST 
dwarf galaxies from both new and comparable quality 
archival imaging. The galaxies range in from —8.23 
(KK 230) to -17.77 (NGC 55) and in distance from 1.3 
Mpc (Sex A) to 4.6 Mpc (DDO 165). We have included 
most galaxies considered dwarf galaxies in the litera- 
ture in our analysis, although the upper mass/luminosity 
cutoff for what con stitutes a dwarf galaxy is some- 
what ambiguous (e.g..lHodgdll97l l: lSkillmanlll996tlMateol 



Il998t ITolstov et al.ii2009D . For example. lMateol (|1998D 
chose to exclude the LMC {Mb = -17.93) and SMC 
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Figure 1. Distribution of the ANGST sample of dwarf galaxies in distance and Ms . Galaxies are color-coded by morphological typo: dSphs 
(rod; T < 0), dis ( blue; T = 10 ), dSpirals (green; T = 8,9), dTrans (cyan; T = 10), and dTidals (magenta; T = 10) (de Vaucoulcurs ct aT] 
ri991l ; iKarachentsev et al.ll2004|) . dTidals and dTrans have been identified b ased on previous literature analysis (see ^2.11 1. We note there 
is likely some ambiguity in the classification of bright dIs and dSpirals (see i|2.1ll . 



{Mb =—16.35) from the 'dwarfs' category. Analo- 
gous to'Mate^ (|1998[) . we have not included NGC 3077 
{Mb =-17.44), NGC 2976 {Mb = -16.77), and 
NGC 300 {Mb = -17.66), which would be among the 
brightest and most massive galaxies in t he sample. De- 
tailed studies of the SFHs of NGC 2976 (iWilliams et al.l 
[2OT0t ) and NGC 300 (jGogarten et al.ll2010[) are available 
in the literature. A full list of sample galaxies and their 
properties are listed in Tabic [T] with the distances and 
blue luminosities of the sample shown in Figure [1] 

Although the sample of ANGST dwarf galaxies is ex- 
tensive, it is not complete within a fixed distance limit. 
Galaxies at low galactic latitudes have been intention- 
ally excluded from the original volume selection to avoid 
complications associated with high degrees of redden- 
ing. While SFHs can still be derived from CMDs o f 
such galaxies (e.g., IC 4662: lMcQuinn et al.|[20Q9ll2Q10l ). 
the effects of extreme reddening can lead to larger un- 
certainties and can require special analysis techniques 
(e.g., individual stellar line of sight reddening correc- 
tions), which detracts from a uniform approach to the 
data reduction. In addition there have been a number 
of recently discovered d warf galaxies in the M81 Group 
(jChiboucas et al.l I2OO90 . which were not discovered in 
time to be included in the original ANGST sample. 



The recent distance reassignment of UGC 4879 to the 
perip hery of the LG (Kopylov ct al. 200^ iJacobs et al.l 
1201 0!) meant this galaxy has also been excluded from the 
ANGST sample. Further, the ANGST sample is likely 
to be missing any faint dSphs which might be located in 
close proximity to MSI, i.e., analogs to Milky Way satel- 
lites such as Draco and Ursa Minor, as such galaxies have 
not yet been detected due to their inherent faintness. At 
a distance of M81 (3.6 Mpc) Draco and Ursa Minor would 
have apparent blue luminosities of 18.98 and 20.13, re- 
spectively. In comparison, Sc22 (m^ = 17.72) has the 
faintest apparent magnitude in the ANGST sample. 
In this paper, we initially divide the sample of 

dwarf galaxies according to morphological typ e, T, 

(|de Vaucouleurs et al.l 119911 : iKarachentsev et al.l 120041 ) 
resulting in 12 dwarf galaxies with T < (dSphs), 5 with 
T = 8 or 9 (dwarf spirals; dSpirals), and 43 with T = 10 
(dIs). For ease of direct comparison with well studied LG 
dwarf gal axies, we hav e adopted nomenclature consis- 
tent with Mateo (1998), in which dwarf ellipiticals (e.g., 
NGC 147, NGC 185) are found to be rare, and dSphs 
are more common. Among the dIs, there could be some 
ambiguity between bright dIs and dSpirals. For consis- 
tency, we defer to the T-type morphological classification 
scheme, but note that the distinction between a bright 
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dl and dSpiral is not always clear. Morphological type T 
= 10 further includes transition dwarf galaxies (dTrans), 
galaxies with reduced recent star formation but high gas 
fractions (e.g., Mated, 1998), and tidal dwarf galaxy can- 
didates (dTidals), which appear to be condensing out of 
tidally disturbed gas. We classify subtypes, i.e., dTi- 
dal and dTrans, as follows: the three dTidals are Holm- 
berg IX, A0952-h069, and BK3N, and are ah located in 
th e M81 Group . For dTrans, we adopt the definition 
of IMateol ()1998l ). namely, that a galaxy has detectable 
gas but very little or no Ha flux. The final sample of 
dTrans was classi fied based on HI and Ha measurements 
in the literature (iCote et allll997HSkillman et"alll2003a : 



Karachentsey et al.l 12004 : Karachentsev fc Kaisin[" 2007 
Begum et all l200a iBouchard et all 120091: iCote et al l 
2009D. We find 12 ANGST dwarf galaxies that satisfy 
the dTrans criteria: KK 230, Antlia, KKR 25, KKH 98, 
KDG 73, ESO294-10, ESO540-30, ESO540-32, KDG 52, 
ESO410-005, DDO 6, and UGCA 438, leaving the final 
tally of true dis at 28. 

Throughout this paper, we adopt the tidal index, 9 
(|Karachentsev et al.l I20Q4) . as a measure of a galaxy's 
isolation. 6 describes the local mass density around 
galaxy i as: 

Q, = max[log{MK/Df^^)] + C, i^l,2,---,N (1) 

where Mk is the total mass of any neighboring galaxy 
separated from galaxy z by a distance of Dik- The val- 
ues of Q we use in this paper have been taken from 
iKarachentsev et all ()2004D . Negative values correspond 
to more isolated galaxies, and positive values represent 
typical group members (see Table [T]). 

2.2. Observations and Photometry 

HST observations of new ANGST targets were carried 
out in two phases due to the failure of ACS in 2007. Prior 
to the failure, we observed new targets using ACS with 
WFPC2 in parallel mode. Galaxies observed post-ACS 
failure were imaged with WFPC2 alone, as part of a 'sup- 
plemental' HST program (Proposal ID s 11307 and 11986 
in Table [n lK~Gilbert et al.llin PreD.h . New ACS obser- 
vations used the F475W (SDSS g') and F814W (/) filter 
combination, to optimize both photometric depth and 
temperature (color) baseline. A third filter in F606W 
(wide V) is also available for most galaxies. The low 
throughput in the bluer filters of WFPC2 led us to only 
use F606W and F814W for WFPC2 obse rvations. 

As described in lDalcanton et all (|2OO90 . both new and 
archival observations were processed uniformly beginning 
with image reduction via the standard HST pipeline. Us- 
ing the ANGST data reduction pipeline, we performed 
photo metry on each image with HSTPHOTQ (Dolphin 
[2000I) . designed for WFPC2 images, and DOLPHOlO, 
running in its ACS-optimized mode, for ACS observa- 
tions, providing for a uniform treatment of all data. The 
resultant photometry for each data set was filtered to en- 
sure the final photometric catalogs excluded non-stellar 
objects such as cosmic rays, hot pixels, and extended 
sources. For the purposes of this paper, we considered 
a star well measured if it met the following criteria: a 

^ http:/ /purcell. as. arizona.edu/hstphot 
^ http:/ /purcell. as. arizona.edu/dolphot 



signal-to-noise ratio > 4 in both filters, a sharpness value 
such that {sharpi + sharp2)^ < 0.075, and a crowd- 
ing parameter such that {crowdi -\- crowd2) < 0.1. To 
characterize observational uncertainties, we performed 
500,000 artificial star tests on each image. Both the full 
and filtered (i.e., the 'gst' files) photometric catalogs, 
HST reference images, and CMDs are publicly available 
on MASTQ. Definitions and detailed descriptions of the 
filtering criteria and the o bservational strategies c an be 
found in Dolphin" (200Cj^ iDaicanton et al.l p009l ). and 
IK. Gil bert et al. (in Prep])- 

Because of the variety of distances in the sample, the 
ACS/WFPC2 field of view does not subtend the same 
physical area for each galaxy. For certain comparisons 
among the sample (e.g., integrated stellar masses), it is 
important to account for these differences in coverage 
area. To that effect, we employ a simple areal normal- 
ization factor based on each galaxy's apparent blue sur- 
face bri ghtness. From the measur ements listed in Ta- 
ble 4 of IKarachentsev et al.l (|2004D . we consider an ef- 
fective elliptical area computed using the angular diam- 
eter and angular axial ratios at a blue surface bright- 
ness level of ^ 25 mag arcsec"^ (^ 26.5 mag arcsec"^ 
in the case of the faintest galaxies). We then calcu- 
late the normalization factor, ^25, for each galaxy by 
taking the ratio of the angular area subtended by the 
ACS / WFPC2 field of view to th e angular area computed 
from IKarachentsev et al.l ()2004D . This normalization has 
been specifically applied to the integrated stellar masses 
throughout this study. The A25 normalizations are typ- 
ically > 1 (i. e., the HST field of view exceeds the area 
computed by IKarachentsev et al.ll200l : see Table |T]), and 
the most and least covered galaxies are BK3N (18.0) and 
NGC 55 (0.02), respectively. 

In most cases, a single HST field was sufficient to cover 
the main optical body of a galaxy, ensuring the SFHs 
are representative of the whole galaxy. However, sev- 
eral of the sample galaxies required multiple observa- 
tions to cover a reasonable fraction of the optical body 
(NGC 2366, Holmberg II, IC 2574, NGC 55, BK5N, 
F8D1, NGC 3109). To derive the SFHs for each of these 
galaxies, we first checked to see that the 50% complete- 
ness limits for each of the fields were similar (i.e., within 
±0.2 mag). This condition was met for all galaxies ex- 
cept IC 2574, NGC 55, and NGC 3109, which have large 
angular sizes and wide ranges of surface brightnesses 
within each galaxy. For galaxies with similar complete- 
ness limits, the photometry and false stars were first com- 
bined, and then the SFH code was run on the combined 
data. For galaxies with overlapping fields (NGC 2366, 
Holmberg II) we carefully removed duplicate stars before 
merging the photometry. IC 2574 was a special case as 
it has three overlapping fields with significantly different 
completeness limits. In this instance, we first removed 
the duplicate stars from overlapping fields, ran the SFHs 
on each field, and then combined the results. 

NGC 55 and NGC 3109 have many new and archival 
HST observations taken over multiple visits. In each 
case, we selected representative non-overlapping fields, 
one in the center of the galaxy and one in the disk. The 
completeness functions were sufficiently different such 
that we did not combine photometry prior to comput- 

"* http:/ /archive. stsci.edu/prepds/angst/ 
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Table 1 

ANGST Dwarf Galaxy Sample 



Galaxy 


Alternate 


Main 


IVlB 


n 


/i-V 


L 




•A2b 


£ lirerh 


^V^FSIAW OU/0 


no 1 


Name 


Name 


Disturber 
















Completeness 


Proposal ID 










(,&} 




I7\ 
\l) 




(.9) 


(^ 

(.10) 


(.11; 


(12j 


KK230 


KKR3 


M31 


-8.49 


1.3 


0.04 


10 


-1.0 


12.0 


F606W,F814W 


-1-0.62 


9771 


BK3N 




M81 


-9.23 


4.0 


0.25 


10 


1.0 


18.0 


F475W,F814W 


-0.32 


10915 


Antlia 




N3109 


-9.38 


1.3 


0.24 


10 


-0.1 


1.20 


F606W,F814W 


+1.75 


10210 


KKR25 




M31 


-9.94 


1.9 


0.03 


10 


-0.7 


2.50 


F606W,F814W 


-0.18 


11986 


FMl 


F6D1 


M82 


-10.16 


3.4 


0.24 


-3 


1.8 


5.00 


F606W,F814W 


+0.11 


9884 


KKH86 




M31 


-10.19 


2.6 


0.08 


10 


-1.5 


5.20 


F606W,F814W 


-1.3 


11986 


KKH98 




M31 


-10.29 


2.5 


0.39 


10 


-0.7 


5.40 


F475W,F814W 


+0.46 


10915 


BK5N-CENTRAL 




N3077 


-10.37 


3.8 


0.20 


-3 


2.4 


3.80 


F606W,F814W 


-0.77 


6964 


BK5N-OUTER 
















3.80 


F606W,F814W 


-0.70 


5898 


Sc22 


Sc-dEl 


N253 


-10.39 


4.2 


0.05 


-3 


0.9 


5.70 


F606W,F814W 


+0.30 


10503 


KDG73 




M81 


-10.75 


3.7 


0.06 


10 


1.3 


12.0 


F475W,F814W 


-0.38 


10915 


IKN 




M81 


-10.84 


3.7 


0.18 


-3 


2.7 


0.58 


F606W,F814W 


-0.82 


9771 


E294-010 




N55 


-10.86 


1.9 


0.02 


-3 


1.0 


4.60 


F606W,F814W 


+1.80 


10503 


A0952+69 




N3077 


-11.16 


3.9 


0.26 


10 


1.9 


1.20 


F475W,F814W 


-0.23 


10915 


E540-032 




N253 


-11.22 


3.4 


0.06 


-3 


0.6 


2.30 


F606W,F814W 


+0.10 


10503 


KKH37 




1342 


-11.26 


3.4 


0.23 


10 


-0.3 


3.70 


F475W,F814W 


+0.12 


10915 


KDG2 


E540-030 


N253 


-11.29 


3.4 


0.07 


-1 


0.4 


2.70 


F606W,F814W 


+0.32 


10503 


UA292 


CVnI-dwA 


N4214 


-11.36 


3.1 


0.05 


10 


-0.4 


5.10 


F475W,F814W 


-0.37 


10915 


KDG52 


M81-Dwarf-A 


M81 


-11.37 


3.5 


0.06 


10 


0.7 


2.30 


F555W,F814W 


+0.37 


10605 


KK77 


F12D1 


M81 


-11.42 


3.5 


0.44 


-3 


2.0 


0.83 


F606W,F814W 


+0.28 


9884 


E410-005 




N55 


-11.49 
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Note. — Properties 


of the sample of ANGST dwarf 


galaxies 


- (1) & 


(2) Names; 


(3) Most 


gravitationally influential 


neighbor (^araehentscv et al.l 



120041) : (4) Absolute Blue Magn itude ( Karae lientsev et al..,2004l'): (5) TRGB Distanee ifDalcanton et al.ll2Mg|: IK. Gilbert et ahilin Prepii): (6) Foreground 



extin ction jSchlegel et al.lll993) : (7) Morphologieal T-Type jde Vaucouleurs et "aTll 19911 : IKarachentsev et al.ll2004f) : (8) Tidal Index I IKaraehentsev et al] 
120041) : (9) Optical Coverage Fraction: Angular area of ACS/WFPC2 coverage divided by angular area of the galaxy measured at a level of ~ 25 mag 
arcsec"^ or ~ 26.5 mag arcsec"^ for 'KK' galaxies ( IKarachentsev et al.|[2004l) : (10) Filter Combination; (11) 50% completeness limit (Mf814w); (12) 
HST Proposal ID, note new ANGST ACS (ID 10915) and WFC2 (11307, 11986) observations. 
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ing the SFH. Instead, we derived the SFHs for each field 
separately, and combined the resultant SFHs from each 
field. The fractional coverage listed in Table [T] takes into 
account the combined areas for galaxies with multiple 
fields. 

3. METHOD OF MEASURING SFHS 

To ensure uniformity in the SFHs o f the ANGST dwarf 
galaxies, we selected one SFH code (jDoTphin 200 2|) and 
one set of stellar evolution models (Padova; Marig o et al.l 
[2008) . The SFH code of Dolphin (2002) provides the user 
with robust controls over critical fixed input variables, 
e.g., IMF, binary fraction, time and CMD bin sizes, as 
well as the ability to search for the combination of metal- 
licity, distance, and extinction values that produce a 
mode l CMD that b e st fits the observed CMD. The mod- 
els of iMarigo et al.l (|2008 h combine updated asymptotic 
gi ant branch (AGE) e volution tracks wit h the models 
of iBertelh et al.l (|1994D (M > 7 M©) and iGirardi et al.l 
(|2002D (M <7 M©). As with any models, there may be 
inevitable systematic biases associated with these par- 
ticular choices (e.g., at some metallicities the Padova 
model RGBs have red offsets from observations; e.g., 
iGallart et al.l 120051 ). However, these systematic effects 
will be shared by all galaxies in the sample, making for 
a robust relative comparison within the sample. We dis- 
cuss systematic uncertainties due to choice of isochrones 
in Appendix IbI 

Here, we briefly summarize the technique of measur- 
ing a SF H based on the full methodology described in 
iDolphinI (11)02). The user specifies an assumed IMF and 
binary fraction, and allowable ranges in age, metallicity, 
distance, and extinction. Photometric errors and com- 
pleteness are characterized by artificial star tests. From 
these inputs, many synthetic CMDs are generated to 
span the desired age and metallicity range. For this work, 
we have used synthetic CMDs sampling stars with age 
and metallicity spreads of 0.1 and 0.1 dex, respectively. 
These individual synthetic CMDs are then linearly com- 
bined along with a model foreground CMD to produce 
a composite synthetic CMD. The linear weights on the 
individual CMDs are adjusted to obtain the best fit as 
measured by a Poisson maximum likelihood statistic; the 
weights corresponding to the best fit are the most prob- 
able SFH. This process can be repeated at a variety of 
distance and extinction values to solve for these param- 
eters as well. 

Monte Carlo tests were used to estimate uncertainties 
due to both random and systematic sources. For each 
Monte Carlo run, a Poisson random noise generator was 
used to create a random sampling of the best-fit CMD. 
This CMD was then processed identically to the origi- 
nal solution, with additive errors in M{,oi and \og(Teff) 
introduced when generating the model CMDs for these 
solutions. Single shifts in Mtoi and log(Te//)) are used 
for each Monte Carlo draw, and the errors themselves 
were drawn from normal distributions with 1-a values 
of 0.41 (Mboi) and 0.03 (log(re//)). These distributions 
were designed to mimic the scatter in SFH uncertain- 
ties obtained by using multiple isochrone sets to fit the 
data, though we prefer to use this technique due to the 
small number of isochrone sets that fully cover the range 
of ages, metallicities, and evolutionary states required to 
adequately model our CMDs. We found that the uncer- 



tainties were stable after 50 Monte Carlo tests, and thus 
conducted 50 realizations for each galaxy. This technique 
of estimating error on SF Hs is described in greater detail 
in I A. DolphinI (|in Prep.l ). 

To minimize systematic effects for comparisons among 
the sample, we selected consistent parameters for mea- 
suring SFHs of all galaxies. All SFHs were measured 
using a single slope power law IMF with a spectral in- 
dex of —1.30 over a mass range of 0.1 to 120 M0, a 
binary fraction of 0.35 with a flat secondary mass distri- 
bution, 71 equally spaced logarithmic time bins rang- 
ing from logt — 6.6-10.15, color and magnitude bins 
of 0.05 and 0.1 mag, and the Padova stellar evolution 
models ([Marigo et al.ll2008l) . We note that the d ifference 
between our selected IMF and a Kroupa IMF ([Kroupal 
1200 It ) is negligible, as the ANGST CMDs are hmited to 
stellar masses > 0.8 Mq. We designated the faint photo- 
metric limit to be equal to the 50% completeness limit in 
each fllter (see Tabled]) as determined by ~ 500,000 ar- 
tiflcial star tests run for each galaxy. The SFH program 
was initially allowed to search for the best fit distance 
and extinction values without constraints. Initial values 
for the dista nces were taken f r om t he TRGB distances 
measured in iDalcanton et al.l (|2009( ). while foreground 
extinction values we re taken from the galactic maps of 
ISchlegel et aT\ (fToOSh . We found no significant discrep- 
ancies between the assumed values and the best CMD fit 
distance and foreground extinction values, i.e., all were 
consistent within error. Final solutions were computed 
usin g the TRGB dist a nces from .Dalcan ton et al.l ()2009f ) 
and IK. Gilbert et al.l din Prep.l ) and foreground extinc- 
tion A^^IueT^^omlSchlere etall (flOQl . 

We placed an additional constraint on the CMD fit- 
ting process, namely that the mean metallicity in each 
time bin must monotonically increase toward the present. 
The deepest ANGST CMDs do not reach the ancient MS 
turnoff, a requisite feature for completely breaking the 
age-metallicity degenera cy of the RGB only using broad- 
band photometry (e.g., iCole et all 120051 : IGallart et al.l 
12005( 1. As a result, SFHs derived from shallow CMDs 
without a metallicity constraint can often have accom- 
panying chemical evolution models that are unphysical 
(e.g., a drop in metalhcity of several tenths of a dex 
over sub-Gyr time scales). A more robust analysis of 
the chemical enrichment of dwarf galaxies would need to 
include the ancient MS, measured gas p hase ab undances, 
and /or individual st ell ar spectra (e . g.. lCole~et al. 200^; 
iMonelli et all l20T0al lbl: IKirbv et al.l 12010^ . We there- 
fore consider analysis of the metallicity evolution of the 
ANGST sample beyond the scope of this paper. 

As an example of a typical CMD flt made by the SFH 
code, we compare the observed and synthetic CMDs of 
a representative ANGST dTrans, DDO 6 (Figure [2]). 
IDalcanton et al.l (|2009f ) determined the TRGB distance 
of DDO 6 to be 3 .31±0.06 Mpc while the foreground 
extinction maps of ISchlegel et al.l (|1998D give values of 
Ab = 0.07 and Ay = 0.06. Allowing the SFH code to 
search for the best fit CMD, we find best fit values of D = 
3.31±0.07 Mpc and Afh^w = 0.05±0.04, both in excel- 
lent agreement with the independently measured values. 
Examining the residual significance CMD, i.e., the dif- 
ference between the data and the model weighted by the 
variance (Panel (d) of Figure [2]) , we see a good, although 
not perfect fit. Notably, the area between the blue helium 



7 



{□) observed CMD (b) modeled CMD 
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Figure 2. A comparison of the observed and model CMDs for ANGST sample dl DDO 6. The observed CMD is shown in panel (a) and 
the best matched model CMD in panel (b). The lower two panels are diagnostic CMDs used to determine the fit quality. Panel (c) shows 
the residual of data — model CMD, with black and white points representing ± 5<j. Panel (d) is the residual significance CMD, which 
is the data— model weighted by the variance in each CMD bin (Dolphin 2002). Based on the residual significance CMD, we see that the 
overall fit is quite good, with only minor discrepancies, most notably between the main sequence and the blue helium burning stars, which 
is likely due to differential extinction effects, but only affects the SFH on time scales < f Gyr. 



burning stars and young MS appears to be too cleanly 
separated in the model, which could be due to differential 
extinction affecting young stars in the observed CMD. 
Additionally, the model red helium burning stars are too 
blue compared to the data, li kely due to uncertai nties in 
the massive star models fe.g.. lGallart et afll2005[ ). How- 
ever, even the most discrepant regions are fit within ± 
5ct, which are indicated by black or white points. Overall, 
the model CMD appears to be in good agreement with 
the observed CM D, ind icating that we have measured a 
reliable SFH. See iDolph in (2002) for a full discussion of 
the quality measures of this CMD fitting technique. 

In general, uncertainties in the absolute SFHs gener- 
ally are somewhat anti-correlated between adjacent time 
bins, such that if the SFR is overestimated in one bin, 
it is underestimated in the adjacent bin. However, cu- 
mulative SFHs, i.e., the stellar mass formed during or 
previous to each time bin normalized to the integrated 
final stellar mass, do not share this property, and thus 
present a more robust way of analyzing SFHs. 

In the following sections, we consider SFHs plotted ver- 
sus two different time binning schemes. For the individ- 



ual galaxies (Figure [3]), we present the cumulative SFHs 
at the highest time resolution possible. Errors in the cu- 
mulative SFHs presented at this time resolution indicate 
both the uncertainty in the fraction of total stellar mass 
formed prior to a given time, and the inherent time res- 
olution for the SFH of that particular galaxy. A coarser 
time binning scheme is used for the absolute SFHs of the 
individual galaxies as well as the mean absolute and cu- 
mulative SFHs across the ANGST sample. In Appendix 
[Xl we determine that six broad time bins of of 0-1, 1-2, 
2-3, 3-6, 6-10, and 10-14 Gyr provide an optimal balance 
between photometric depths and age leverage contained 
in the CMDs of the ANGST dwarfs. The broader time 
bins allow us to securely average the best fit SFHs from 
individual galaxies. For the absolute SFHs of individual 
galaxies, broader time bins tend to minimize co- variance 
between adjacent time bins. Thus applying a broader 
time binning scheme to the absolute SFHs of individual 
galaxies provides for more secure measurements of the 
SFRs. 
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Figure 3. The SFHs and cumulative SFHs, i.e., fraction of total stellar mass formed prior to or during a given time bin, of 
ANGST sample galaxies presented in order of increasing absolute blue luminosity. The grey dashed line in the SFH plots is the 
lifetime averaged SFR. In the cumulative SFHs, this same rate is also represented by the grey dashed line with a slope of unity, 
i.e., a constant SFH. The error bars (or yellow envelope for the cumulative SFHs) shown are for the 16th and 84th percentile for 
the distribution of SFHs as computed via the Monte Carlo process described in iJ21 The axes of the SFHs have been scaled so 
that galaxies of comparable luminosity are on similar scales. Some galaxies absolute SFHs have been scaled up to clarify details. 
Note that some of the SFRs have been scaled for clarity as indicated in the absolute SFH panel (e.g., '5x' means the SFH has 
been multiplied by a factor of 5 or '0.25x' means it has been scaled down by a factor of 4). The optical coverage fraction (see 
Table [T]) indicates how much of the optical galaxy a SFH represents. 
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Figure 3. Continued 
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Figure 3. Continued 
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Figure 3. Continued 
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Figure 3. Continued 
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Figure 3. Continued 
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Table 2 

SF Properties of the ANGST Dwarf Galaxies 
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Note. — SF properties of the ANGST dwarf galaxies from the CMD-based SFHs - (1) Galaxy Name; (2) Lifetime 
averaged SFR, i.e., average SFR over the history of the galaxy, from the best fit SFH; (3) Integrated stellar mass from 
the best fit SFH, with the A.ih normalization applied; (4) Indicative mass-to-light ratio in solar units; (5)— (10) The 
cumulative fraction of stars formed prior to 10, 6, 3, 2, 1 Gyr ago from the best fit SFHs. 
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4. THE EXTENDED STAR FORMATION HISTORIES OF 
DWARF GALAXIES 

The ANGST dwarf galaxies exhibit a wide variety of 
complex SFHs. In Figure [3] we show the absolute SFHs, 
i.e., SFR(t), and the cumulative SFHs of the individ- 
ual galaxies, sorted in order of increasing blue luminos- 
ity. A cursory inspection of the 60 SFHs reveals that 
often galaxies with similar luminosities, morphologies, 
or chemical compositions, do not have consistent SFHs, 
confirming the complexity of dwarf galaxy SFH s pre- 
viously found in studies of the LG (e .g., iMated [19981 : 
iDolphin et all [20051: l^stov et al.|[200a . The wide vari- 
ety of SFHs underlines the importance of having a large 
sample; conclusions based on the SFH of a single galaxy 
may not be be representative of the population as a 
whole. 

This paper leverages the size of the ANGST sample to 
explore general trends seen in dwarf galaxies, e.g., differ- 
ences among the morphological types. While the SFHs 
of individual galaxies can provide useful insight into a 
particular system's evolutionary state, we generally do 
not consider individual galaxies in this paper. However, 
in Table [2] we provide data on the SFHs of individual 
galaxies considered in this sample, which may be useful 
for specific studies. 
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Figure 4. The mean specific SFHs (sSFH), i.e., the SFH divided 
by the integrated stellar mass, for each morphological type. The 
error bars reflect the uncertainties in the mean sSFHs. Generally, 
dwarf galaxies had a high level of SF at ancient times (> 10 Gyr 
ago), with a low level of SF at intermediate times (1-10 Gyr ago). 
There are noticeable differences in the sSFHs among the morpho- 
logical types within the most recent 1 Gyr. 



We first consider the unweighted mean specific SFHs 
(i.e., a SFH divided by the total integrated stellar mass 
formed) of the ANGST sample grouped by morphologi- 
cal type (FigureH]). Comparison among the specific SFHs 
shows general consistency for times > 1 Gyr ago, with 
the exception of dTidals (note that dTidals will be dis- 
cussed in greater detail i n the context of recent SFHs in 
ID. Weisz et al..|[in Prep.l ). Qualitatively, we find that a 
typical dwarf galaxy exhibits dominant ancient SF (> 10 
Gyr ago), and lower levels of SF at intermediate times 
(1-10 Gyr ago). The only consistent difference among 
the morphological types is within the last 1 Gyr, where 
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Figure 5. The individual and mean cumulative SFHs, per mor- 
phological type. Error bars have been omitted for clarity, and are 
instead listed in Tables [2] and |3] The dot-dashed line represents 
50% of the total stellar mass. While the individual cumulative 
SFHs show a wide variety, the mean cumulative SFHs are remark- 
ably similar. Most dwarf galaxies are not entirely old stellar pop- 
ulations, i.e., they have intermediate or recent SF, and none are 
consistent with simple models of SF, i.e., single epoch or constant 
SFHs. Excluding dTidals, the average dwarf galaxy formed > 50% 
of its stars prior to z ~ 2 (10 Gyr ago). 



dSphs exhibit a significant drop in SFRs relative to the 
other types. This suggests that many of the dSphs in 
the ANGST sample could have been gas-rich as recently 
as 1 Gyr ago, implying that the process of gas loss can 
occur relatively quickly and at late times. It is tempting 
to associate the sharp drop in the SFR of a typical dSph 
~ 1 Gyr ago with rapid gas loss. However, broad un- 
certainties in the A GE models (e.g.. iGirardi et al.l[2010l : 
iMarigo et al]|2010[) make the precise age for the drop in 
dSph SFRs and the degree of synchronization uncertain. 
In addition, subsamples of the dSphs do appear to have 
dramatically lower SFRs at more intermediate ages. 

Cumulative SFHs allow us to readily compare galaxies 
of different masses. Like the absolute SFHs, the cumula- 
tive SFHs plotted in Figure [5] show significant diversity 
within each morphological class, yet converge on mean 
values that are broadly consistent among the morpho- 
logical types. Perhaps the most striking result is that 
the average dwarf galaxy formed > 50% of its total stel- 
lar mass by z 2 (10 Gyr ago) and > 60% by z ~ 1 (7.6 
Gyr ago), independent of morphological type (Table (H]). 

Although dSphs appear to have formed a slightly larger 
percentage of total stellar mass than dis by z ~ 2, 
additional consideration of systematic uncertainties, as 
discussed in Appendix |b1 indicate that the ANGST 
data is not deep enough to make such fine distinctions. 
More importantly, the amplitude of the difference is sig- 
nificantly smaller than typical models, which assume 
that dSphs are dominated by ancient stellar popula- 
tions and dIs are consistent with constant SF over their 
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Table 3 

Mean Star Formation Properties per Morphological Type 



Group 


(M25> 

(107 M0) 


flOGyr 


fdGyr 


hOyr 


hOyr 


hOyr 


fCnrreut 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



dSph 
dl 

dTrans 
dSpiral 
dTidal 



5.4±0.9 
7.3±1.8 
2.8±1.5 
74±35 
1.6±0.35 



0.65 
0.54 
0.53 
0.59 
0.38 



-0.06 
-0.08 
+0.05 
0.04 
+0.10 
0.09 
+0.14 
0.11 
+0.08 
0.05 



0.71 
0.67 
0.63 
0.71 
0.64 



+0.05 
0.06 
+0.03 
0.04 
+0.09 
0.08 
+0.10 
0.06 
+0.09 
0.05 



0.87tH^ 
0.86lH^ 

+0.03 
-0.02 
,-, OO+0.06 



0, 



0.89 
0.82 
0.89 
0.89 
0.84 



+0.03 
0.02 
+0.02 
0.01 
+0.01 
0.02 
+0.03 
0.02 
+0.07 
0.07 



0.98 
0.92 
0.96 
0.95 
0.85 



+0.00 

0.00 
+0.01 

0.01 
+0.01 

0.01 
+0.01 

0.01 
+0.08 

0.06 



1.00 
1.00 
1.00 
1.00 
1.00 



Note. — Mean SF properties for the ANGST dwarf galaxies, grouped by morphological 
type. (2) The integrated stellar mass normalized to the A25 areal fraction and then averaged 
per morphological type; (3)— (7) The cumulative fraction of total stellar mass formed prior to 
10, 6, 3, 2, 1 Gyr ago. Error bars represent the uncertainty in the mean value. 



lifetimes (e.g lBinggelil 119941: ISkillman fc Bended IT995I: 
IGrebel et al.l[2003l) . 

The mean cumulative SFHs hint at divergence among 
the morphological types within the last few Gyr (Figure 
[5]). However, the exact characteristics of differences in 
the SFHs during this time period are ambi guous. Broad 
uncertainties in the AGB star models (e.g.. lGirardi et aD 
I2OIOI: iMarieo et all l2010l ). our modest time resolution, 
and systemic effects on the mean cumulative SFH (see 
Appendix |B| all affect our ability to quantitatively ex- 
amine the details of intermediate age SF. 

We do identify secure measurements of differences 
within the most recent 1 Gyr. At these times, lumi- 
nous MS and core blue and red h elium burning stars 
provi de excellent age leverage (e.g.. iDohm-Palmer et al.l 
Il998f) . and the SFHs clearly illustrate differences among 
the morphological types. Specifically, the typical dSph, 
dl, dTrans, and dSpiral formed ~ 2%, 8%, 4%, and 5% of 
their total stellar mass within the most recent 1 Gyr. The 
combined findings from intermediate and recent SFHs 
suggest that morphological differences and complete gas 
loss in dwarf galaxies can be relatively recent phenomena, 
at least within the Local Volume. dSphs that are strictly 
old, i.e., have no AGB star populations, are known to 
exist in the LG (e . g.. iMould et al1ll982t IMateol 119981: 
IDolphin et all 120051: rfolstov et al.ll2009D . but appear to 
only represent a minority of all known dSphs in the larger 
volume we analyze here. This is at least in part due to 
the detection limit of the ANGST sample as discussed in 

To further explore the differences among the morpho- 
logical types, we consider the fractional SFH, i.e., the 
fraction of stellar mass formed in each time bin, as a 
function of the total normalized stellar mass, M25 (com- 
puted by integrating the SFH over time and applying the 
A25 area normalization). We show values for both the 
individual galaxies (grey symbols) and the mean values 
per morphological type, with error bars representative 
of the uncertainty in the mean in Figure [S] The mean 
stellar masses and cumulative SFHs for each of the mor- 
phological types are listed in Table |3l 

The mean stellar masses show that dSpirals are typ- 
ically the most massive galaxies, dTrans are the least 
massive, and dis and dSphs have similar mean masses. 
There are not strong links between patterns of SF, total 
stellar mass, and morphological type, confirming that SF 
processes among the different types of dwarf galaxies are 



not dramatically different. We continue analysis involv- 
ing the integrated stellar masses in ij5.ll 

5. THE MORPHOLOGY-DENSITY RELATIONSHIP 

Galaxies in clusters, groups, and the field follow simi- 
lar morphology-density relationships. Namely, gas-poor 
galaxies, i.e., ellipticals, are generally found t o be less 
isolated than gas-rich galaxies, i.e., spira ls (e.g.. iDressled 
[19801 IOemlei<ll97l IBlanton et aDf2005l) . 

This same morphology-density r elationship ha s been 
found for dwarf galaxies (e.g., 'Einas to et al.l 11974 : 
Mateo 1998; van den Bcrgh 2000; Skillm an et 111120033 : 
Geha et al.l 120061 ). The change in the relative fraction 



of gas- rich and gas-poor galaxies with environment, pro- 
vides a simple test for various models of dwarf galaxy 
evolution. As illustrated in Figures [7] and HI the ANGST 
dwarf galaxies clearly adhere to the morphology-density 
relationship, when using tidal index, Q, as a proxy for 
local density. Typically, dIs are significantly more iso- 
lated than dSphs, in spite of having similar mean total 
stellar masses. dTrans, on average, have intermediate 
tidal indices between dIs and dSphs (see ^5.'6\i , and have 
a lower mean stellar mass than either. dSpirals are the 
most massive galaxies, and are typically located in re- 
gions of intermediate isolation. These findings are in 
general a.greement with earlier studies of LG dwarf galax- 
ies_(e.g.. lMateolll998l:lvan den Berg"hll2000l: l^stov et al.l 
I2OO9I) . 

When combined with measured SFHs, the 
morphology-density relationship can be used to 
gauge the reliability of dwarf galaxy evolution models, 
particularly processes that induce gas loss. Scenarios 
favoring internal mechanisms (i.e., stellar feedback) as 
the primary driver of gas loss (e.g.. iDekel &: Sill3 119861 : 
iDekel fc Wool 120031) can reproduce a number of ob- 
served dwarf galaxy properties (e.g., surface brightness , 
rotation velocities, metallicities, etc; iWoo et alTl2008| ). 
However, such models are generally unable to account 
for t he morphology-density relationship (e.g., iMaveii 
|2010[) and often predict that gas-rich and gas-poor dwarf 
galaxies may have different patterns or efficiencies o f 
SF (e.g.. IDekel fcSili3[T986t Skillma n fc Bended [T995h . 
In contrast, some models that additionally factor in 
external effects (e.g., gravitational interactions, ram 
pressure stripping) have been able to reproduce a wide 
range of dwarf galaxy properties, includin g a canonical 
morphology-density re lationsh ip (e - g.. |Maver et all 
I200I iMaveil l20Ta iKazantzidis eTaLl I20T0I) . In the 
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Figure 6. The fractional SFHs, i.e., tiie fraction of total stellar mass formed during a given time bin, for each galaxy shown as a function 
of the total stellar mass and morphological type. M25 is derived by integrating the SFHs (Figure [Sj, and normalizing to A25 to account for 
differences in the observed areas. The time bins are: /o-iG = O-lj /l— 2G = 1-2, /2— 3G = 2-3, fs—eG = 3-6, /e— lOG = 6-10, /lo— 14G = 
10-14 Gyr ago. The grey points represent individual galaxies, while the colored points show the mean values per morphological type with 
error bars representing the uncertainty in the mean. 



following sections, we discuss results from the ANGST 
SFH analysis within the context of physical processes 
that can affect the evolution of dwarf galaxies. 

5.1. Comparing Stellar and Gas Masses 

The relative masses of the gas and stellar components 
can provide clues to dwarf galaxy evolution. For the 
ANGST sample, we consider the gas mass and bary- 
onic gas fraction as functions of total stellar mass and 
ti dal index. The gas masses a re based on the HI masses 
in IKarachentsev et all (j2004[ ). corrected by a factor of 
1.4 to account for helium content, while the baryonic 
gas fractions, Mgas/ Mbaryonic, is defined to be Mgas/ 
{Mgas+M25), and does not account for the contribution 
due to warm/hot baryons. The stellar masses, M25, are 
derived from integrating the SFHs over time and apply- 
ing the areal normalization, A25 ■ 



We first consider the relationship between total gas 
mass and integrated stellar mass as shown in Figure 
[HI where the dot-dashed line denotes Mgas = M25. In 
this context, the most massive galaxies occupy the up- 
per right portion of the plot, while lower mass dwarfs are 
located to the left. Galaxies without detectable gas are 
located in the lower portion of the plot, and have been 
placed at \og{Mgas) = 3 for convenience. To illustrate 
the range of stellar mass covered by the ANGST dSphs 
relative to the LG dSphs, we have included select LG 
dSphs in gr ey (pl aced using stellar mass-to-light values 
from Mateol ll998D . 

Figure|9]provides a concise snapshot of the current evo- 
lutionary state of nearby dwarf galaxies. For large stellar 
masses, this view illustrates the morphological ambigu- 
ity between dis and dSpirals; galaxies from both classes 
are gas-rich and can have stellar masses > 10^ M0. HI 
surface density maps of some dIs even reveal hin ts of spi- 
ral gas structure in dIs (e.g., iPuche et al.lll992l) . In this 
same mass regime, however, there is a conspicuous ab- 
sence of dSphs. In the ANGST sample, we not do find 
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Figure 7. The observed morphology-density relationship for the ANGST dwarf galaxies. Values of M25 have been derived by integrating 
the SFHs in Figure [ 3l and normalizing to ^25 to account for differences in the observed areas. The tidal indices, ©, have been taken from 
[Karachcntse v et al.l 1120041) . Negative values of represent isolated galaxies, while positive values represent typical group members. 



dSphs more massive than ~ 10* M©, which is in agree- 
ment wit h the lack of massive gas-poor g alaxies in the 
LG fe.g- iMateol [19981: PrIsFoy et alj 120091) . dTrans are 
predominantly located at low stellar masses, yet have rel- 
atively large gas supplies. These properties suggest that 
some dTrans may not be significantly different from low 
mass dis; we return to this point in §5.31 

The baryonic gas fraction provides a more detailed 
view of the current evolutionary state of the ANGST 
dwarf galaxies (Figure [TU| . The general trend for dwarf 
galaxies mirrors t hat of massive counterpar ts on the Hub- 
ble sequence fe.g. JRoberts &: Havnesl[T99l ) . namely that 
spiral galaxies typically have low gas fractions, while 
irregulars have higher gas fractions. This perspective 
also reinforces the view that many dTrans could be 
low mass dIs in betw een episodes of massive SF (e.g., 
iSkillman et"alll2003al) . 

There appear to be several outliers to the general 
trends in Figure [TUl The two most conspicuous outliers, 
HS 117 and DDO 113, may simply be morphological mis- 
classifications. HS 117 is classified as a dl, but has a very 
low gas content, with an upper limit of Mhi ~ 10^ M© 



(jHuchtmeier et al.l[2000l) and error bars consistent with 
zero. Interestinglv. iKarachentsev fc Kaisinl ()2007l ) detect 
low levels of Ha in HS 117, which reinforces the dl clas- 
sification, but based on the lack of HI, they deem this 
a dSph. Morphologically, inspection of the HST image 
further reveals that it appears to be superimposed on 
an HII region, leading to an erroneous dl classification. 
We therefore suggest HS 117 is best classified as a dSph 
or possibly a dTrans. DDO 113 is likewise classified as 
a dl with negligible gas content, but ha s a possible Ha 
detection (Karachc ntsev fc Kaisinl [20071). Inspection of 
the HST based CMD suggests that DDO 113 resembles 
a prototypical dSph, with a handful of relatively faint 
blue stars, suggesting it is also likely a dSph or a dTrans. 

Additional putative outliers include three dTrans 
(ESO294-010, ESO410-005, ESO540-032) and three dIs 
(NGC 4163, IC 5152 and DDO 125). Two of the 
dIs (IC 5152 and DDO 125) are isolated and relative 
massive, yet have moderate gas fractions. The three 
dTran s have confirmed blu e horizontal branch popula- 
tions (|Da Costa et al.|[2010f) , which helps constrain their 
ancient SFHs. However, none of these galaxies show any 
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Figure 8. The dwarf galaxy morphology-density relationship represented by the fractional SFH, i.e., the fraction of total stellar mass 
formed during a given time bin, as a function of tidal index, 0, and morphological type. The time bins are such that fo —io = 0-1, .fi—2G 

' f 2—sG = 2-3, fs—eG = 3-6, fe—ioG = 6-10, /lo— 14G = 10-14 Gyr ago. The tidal indices, ©, have been taken from IKarachentsev et ahl 
II2004I ) . Negative values of © represent isolated galaxies, while positive values represent typical group members. The grey points represent 
individual galaxies, while the larger colored points are the mean values. Error bars are the uncertainties in the mean values. 



unique features in their SFHs that would explain their 
outher status. In contrast, NGC 4163 is a starburst dwarf 
galaxy, which may have consumed a significant amount 
of gas du e to an intense burst o f SF within the last ~ 1 
Gyr fe.g- lMcQuinn et al.|[2009l ). 

Examining the gas fraction as a function of isolation 
(Figure fTTj) . we see that isolated galaxies tend to have 
high gas fractions, while those in high density environ- 
ments have low gas fractions. Interestingly, there appear 
to be no galaxies with any appreciable gas in dense envi- 
ronments {& > 1.5; HS 117 and DDO 113 have uncertain 
gas measurements as described above). We also see some 
evidence for separation of dTrans into groups of isolated 
(8 < 0) and moderate/high {Q < 0) density environ- 
ments, which we further discuss in i i5.3l 

5.2. Mechanisms for Complete Gas Loss 

How gas-rich galaxies completely lose their gas has long 
been an outstandi ng question in dwarf galaxy evolution 
(e.g., lMaveii[2010l and references therein). Within the 
ANGST sample, we can test the viability of various gas 



loss mechanisms using the observed properties of gas- 
poor dSphs. The ANGST dSphs share several notable 
characteristics including: (1) little SF in the most recent 
1 Gyr compared to gas-rich dwarfs, (2) similar total stel- 
lar masses, (3) SFHs that are generally extended and in- 
distinguishable from dis, and (4) are located exclusively 
in high density environments. In what follows, we con- 
sider the impact of putative gas loss mechanisms on the 
evolution of a typical gas-rich dwarf galaxy (e.g., with 
both AI25 and Mgas ^ 10^ M©) in the context of Figure 

E 

The first mechanism for gas removal is consumption 
through SF. This process will increase the stellar mass 
and decrease the gas mass of our prototypical dl, moving 
it to the right on Figure El Though SF can consume large 
amounts of a galaxy's gas reservoir, the gas densities will 
eventually becom e too low to continue to form stars (e.g., 
lKenmcutqll989[ ). This suggests that if consumption was 
the only mechanism, we should observe trace amount of 
gas in dSphs, which is not typically the case. 

In the event that gas density was not a limiting fac- 
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Figure 9. Log(iV/Gas) plotted versus Log(Af25), where M25 has been derived by integrating the SFHs (Figure [3l), and normahzing 
\0A25 to account for difTerences in the observed areas. We computed Moas by correcting the HI masses from lKarachentsev et al.l 
1)20041 ) by a factor of 1.4 to account for hehum content. The dot-dashed hue represents the Mcas = M25 equaUty. Galaxies 
without detectable gas are located in the lower portion of the plot, and have been placed at log(Mgas) = 3 for convenience. To 
illustrate the range of stellar masses spanned by the ANGST sample, we have included select LG dSphs in grey. Considering 
a present day gas-rich galaxy, SF will increase the stellar mass and decrease the gas mass, moving the galaxy down and to the 
right. Gas removal will move a galaxy downward, while gas addition, e.g., accretion, moves a galaxy up. Stellar mass loss moves 
a galaxy to the left. SF and stellar feedback alone cannot transform a gas-rich galaxy into a dSph, however models including 
external effects, e.g., ram pressure stripping and tidal forces, can remove sufficient stellar and gas mass to transform a gas-rich 
into a gas-poor galaxy (see i ^5.2p . 



tor, characteristically low SFRs and SF effi ciencies in 
dwarf galaxies (~ 1%; e.g.. lLerov et al.l[2008l ) imply long 
gas consumption timescales. In this event, gas supplies 
would typically not be exhausted for more than a Hub- 
ble time, which suggests that dSphs would be quite rare. 
In addition, gas consumption does not have a depen- 
dence on the environment external to a galaxy. Thus, 
it cannot account for the observed morphology-density 
relationship. Although SF is a critical process to galaxy 
evolution, it cannot be responsible for the complete re- 
moval of gas from dwarf galaxies. 

The second possible mechanism for removing gas is 
stellar feedback. Mechanical energy due to stellar winds 
and supernovae provide an appealing expla nation for gas 
removal from a galaxy via expulsion (e.g.. iDekel fc SiI3 



Il986f ). For the prototypical gas-rich dl, this process re- 
moves gas with minimal effect on stellar mass, moving 
the galaxy straight down in Figure [5] Several feedback 
models can reproduce observed dwarf galaxy relation- 
ships between mass, lum i nosity, and rotational velocities 
(e.g. IDekel fc Silkl 119861 : fPekel fc Wool[200l IWoo et^ 
120081 ). Additionally, these models are generally able 
explain the observed mass-mctallicity relationship for 
dwarf galaxies (e.g., [Lee et al. 2006). 

Although these models hold promise, there are two sig- 
nificant challenges to gas removal in dwarf galaxies due 
to feedback. First, a number of simulations have demon- 
strated that the energy due to stellar feedback is insuf- 
ficient to completely expel cold gas (e.g., HI) from the 
gravitational potential of a typical dwarf galaxy (e.g.. 
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normalized to A25 to account for differences in the observed areas. The dot-dashed line represents the Moas = JV/25 equality. 



Mac Low fc Ferraral fl99a iD'Ercole fc Brighentil 11991 
Marcolini et al. l200a iRevaz et s\\\2Q0 9]~ Second, stellar 
feedback alone cannot explain the morphology-density 
relationship. Similar to the rationale for gas consump- 
tion through SF, stellar feedback is not dependent on 
the environment surrounding a galaxy. Thus, while SF 
and stellar feedback are primary drivers of dwarf galaxy 
evolution, it is unlikely that either or both are able to 
account for the final transition from a gas-rich to a gas- 
poor state. 

Th e third mechanism w e consider is ram pressure strip- 
ping. iLin fc Fabeii ()1983D first proposed that dSphs were 
once gas-rich dwarfs whose gas supply has been removed 
by ram pressure stripping from the intergalactic medium 
(IGM). This mechanism has negligible impact on stel- 
lar mass loss, thus moving our prototypical dl straight 
down in Figure |9l Satellite galaxies in higher density en- 
vironments would likely encounter a denser IGM, which 
would lead to more efficient gas loss. This provides a 
feasible explanation for the observed morphology-density 
relationship. 



Despite the appeal of this mechanism, there is a dis- 
tinct lack of observational evidence of systematic ram 
pressure stripping o f LG dwarf gal axies. For exa mple, 
iMcConnachie et all (2007) and Kniaz ev et all (|2009 ) find 
that ram pressure stripping appears to be a localized 
phenomena, and only evident for a small minority of 
LG dwarfs. However, given that the LG is in a state 
of relatively passive evolution, it could be that signa- 
tures of ram pressure stripping have been erased in many 
satellites. An additional challenge c omes from th e mag - 
nitude of ram pressure stripping. iMaver et al.l ()2006f ) 
demonstrate the ram pressure stripping is unlikely effec- 
tive enough to completely remove the gas supply of a 
dwarf galaxy in the LG. It appears that ram pressure 
stripping is not able to account for complete gas loss in 
dwarf galaxies. 

We next consider tidal effects on gas removal from 
dwarf galaxies. The close passage of a dwarf galaxy to 
a massive companion has strong gravitational effects on 
the gas, stellar, and dark matter conten ts of the smaller 
galaxy. As shown in lMaver et al.l (|2001bf ). the magnitude 
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Figure 11. The ratio of gas mass to total b aryonic mass (i.e., the gas fraction; Moas to Mcas + M25) plotted versus the tidal index, 
©. The tidal indices, ©, have been taken from lKarachentsev et aP ||2004). Negative values of O represent isolated galaxies, while positive 
values represent typ ical group members. HI m asses have been taken from Karachontsev ot al. (2004). We computed Mq^s by correcting 
the HI masses from lKarachen tsev et aLl 1)20041 ') by a factor of 1.4 to account for helium content. The total integrated stellar masses are 
from the SFHs (Figure |3}, normalized to A2Z to account for differences in the observed areas. Gas-poor galaxies typically have positive 
tidal indices, while gas-rich galaxies have predominantly negative values. It is interesting to note that for © > 1.5, there are no gas-rich 
galaxies. 



of the tidal force during an interaction does not appear 
to be enoug h to remove t he gas conte nt from a gas-rich 
dl. Instead, iMaver et al] (|2001aL 12006 ^) advocate a more 
complex approach in which a combination of ram pres- 
sure stripping and stellar mass loss due to tidal effects 
provide a plausible model, 'tidal stirring', for the trans- 
formation of gas-rich dis into gas-poor dSphs. In this 
scenario, the prototypical gas-rich dl is able to lose both 
stellar and gas mass, allowing it to move down and left 
in Figure ini meaning a dl with a high stellar mass could 
be transformed into a less massive dSph. 

Predictions from tidal stirring models provide quali- 
tative explanations for a number of trends seen in the 
ANGST sample. Foremost, tidal stirring naturally pro- 
duces a morphology-density relationship. Galaxies in 
lower density environments have had few (or no) inter- 
actions with massive galaxies, and are able to maintain 
high gas fractions. This is in general agreement with the 
trends seen in Figure [TT] Of the galaxies that do in- 



teract with a massive companion, tidal stirring predicts 
that stellar and gas mass loss happens progressively over 
several Gyr (a single complete orbit in the LG is typi- 
cally 1-2 Gyr; Kazantzidis et al. 2010). Therefore, these 
galaxies spend the majority of their life cycling between 
SF and stellar feedback. By extension, this implies that 
many dSphs and dIs should not have drastically different 
SFHs, in general agreement with the SFHs of ANGST 
dwarf galaxies. 

While the general qualitative agreement between tidal 
stirring predictions and the ANGST sample is encourag- 
ing, a quantitative comparison between predictions and 
observed dwarf galaxy SFHs are needed to place precise 
constraints on the gas loss processes in dwarf galaxies. 

5.3. The Nature of Transition Dwarf Galaxies 

With little evidence of recent SF, yet detectable 
amounts of HI, dTrans may hold clues to the transfor- 
mation of gas-rich to gas-poor galaxies. Synthesizing 
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several previous studies, there are two favored scenar- 
ios for the origins of dTrans, namely that they are ei- 
ther in the last throes of SF or are observed during tem- 
porar y lulls between episodes of massive star forrnation 
(e.g., [ Cote ct al' ^199?: 'Mateo 199^: 'Gre bel et al] [200l 
iSkillman e t al. 2003a; Cote ct al. 2009). 

dTrans in the ANGST sample have characteristics 
consistent with both scenarios. Those with low gas 
fractions (MgaJMbaryomc < 0.4; ESO294-010, ES410- 
005, and ESO540-032) typically have positive tidal in- 
dices, indicating a higher likelihood of gas loss via in- 
teractions (Figure [TT|) . These three galaxies may be 
examples of galaxies genuinely transitioning from gas- 
rich to gas-poor states. In contrast, the most isolated 
dTrans (KK 230, KKH 98, KKR 25, and UA 438) typi- 
cally have high gas fractions (> 0.9) and stellar masses 
among the lowest in the ANGST sample (Figure [T0|) . 
These dTrans have properties consistent with other low 
mass dis, and could simply be low mass dis lacking 
Ha, which can be attributed to a wide range of ef- 
fects including t emporary episodes bet ween massive star 
formation (e.g.. ISkillman et "ani2003af ). a stoc hastically 
sampled or systematically var ying IMF (e.g., lEldridgd 
[Ma iWeidner fc Kroupall2005h. and /or leakage of ioniz- 
ing photons (e.g.. [Ferguson "etall [19961) . The remaining 
dTrans (AntUa, KDG 2, DDO 6, KDG 52, KDG 73) 
appear to be have masses and tidal indices close to the 
genuinely transitioning group, but gas fractions similar 
to the low mass dis. It is possible that these galaxies 
may have had few interactions, resulting in less gas loss. 
Detailed multi- wavelength studies of these galaxies could 
provide insight into gas loss mechanisms in dwarf galax- 
ies. 

6. DWARF GALAXIES IN A COSMOLOGICAL CONTEXT 

As the smallest and most pristine galaxies in the Uni- 
verse, dwarf galaxies occupy a critical, but poorly con- 
strained role in cosmol ogical models of galaxy forma- 
tion and evolution (e.g.. lMoore et al.lll999i : [Klvpin et al.l 
[1999). Because of their intrinsic faintness, dwarf galax- 
ies have been difficult to d irectly detect at cosm ologically 
significant redshifts (e.g.. iMathews et "aril2004[ ). 

As an alternate approach to observing high redshift 
dwarf galaxies, we present a comparison between the 
mean cumulative SFHs of the ANGST sample dwarfs 
and the cosmic SFH, as derived from observed UV fluxes 
in high redshift galaxies (e.g., Reddy et al, 2008^ . In the 
top panel of Figure [121 we see good agreement between 
the cosmic SFH (grey shaded region) and the mean cu- 
mulative SFHs of the ANGST dwarf galaxies at z ^ 2 
(10 Gyr ago). At more recent times, the cosmic SFH 
increases sharply, while the cumulative SFHs of dwarf 
galaxies increase at a slower rate. 

This comparison suggests that the largest differences 
between SFRs in dwarf galaxies and massive galaxies 
studied at higher redshifts are at intermediate and re- 
cent epochs. This finding is generally consistent with 
the expectations of galaxy 'downsizing', where high mass 
galaxies pr eferentially stop fo rming stars at higher red- 
shifts fe.g.. lCowie etalllTggl . The precise role of dwarf 
galaxies in theories of downsizing is highly uncertain, 
making this result diffic ult to interpret in the co ntext 
of current models (e.g., iMouri fc Taniguchj [20061) . We 
also caution that the physical significance of this find- 
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Figure 12. The mean cumulative SFHs for dSphs (red), dis 
(blue), and dTrans (purple). The error bars represent the 
uncertainties in the mean , while the dot-dashed line indi- 
cates 50% of the total stellar mass. In the top panel, the 
Sr ey shaded reg i on re presents the cosmic SFH as measured 
by [Reddy etUT pOOl ). We see general good agreement be- 
tween the lifetime SFHs of dwarf galaxies and the cosmic 
SFH, with some discrepancies at intermediate times (see ^ 
for further discussion). In the bottom panel, the grey shaded 
region represents the expected cumulative SFH for an expo- 
nentially declining SFH model, for a range of r values from 
0.1 to 14.1 Gyr. Notice that single valued r models are not 
good representations of the observed SFHs. 



ing is still ambiguous due to a several factors including 
uncertain AGB star models and redshift d ependent dust 
corrections applied to the cosmic SFH (e.g..lGirardi et al.l 
[MotlMarigo et al.l[20Tol: IReddv et al]l2010D . A compar- 
ison of SFHs from the entire ANGST sample (including 
the most massive galaxies) and the cosmic SFH is pre- 
sented in B. Williams et aP (|in Prepf ) . 

The cumulative SFHs can also provide insight into the 
validity of model SFHs often used to describe the evo- 
lution of dwarf galaxies. As an illustrative example, we 
have selected a simple exponentially declining SF model 
(SFR oc e"*/"^; a r model), and show the cumulative 
SFHs for a range of t values (grey shading) in the lower 
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panel of Figure [T^l Visual inspection suggests that the 
measured SFHs show significant deviations from the pre- 
dicted smooth curve of single r valued models. A test 
between the measured cumulative SFHs and exponen- 
tially declining SFH models (with r varying from 0.1 to 
14.1 Gyr) confirms that any single value of r does not 
accurately represent the data. Similarly, the mean mea- 
sured SFHs also are not well matched with a constant 
SFH (i.e., T — > oo Gyr) or a single ancient epoch of SF 
followed (t < 0.1 Gyr) by passive galaxy evolution. In- 
stead more complex and multi-component models may 
better describe the lifetime SFHs of dwarf galaxies. 

7. SUMMARY 

We uniformly analyzed SFHs of 60 dwarf galaxies in 
the nearby universe based on observations and data pro- 
cessing done as part of the ANGST program. While the 
SFHs of individual galaxies are quite diverse, we find that 
the mean SFHs of the different morphological types are 
generally indistinguishable outside the most recent ~ 1 
Gyr. On average, the typical dwarf galaxy formed ~ 50% 
of its stellar mass by z '--^ 2 and 60% by z ~ 1. 

Among the morphological types, the SFHs hint at di- 
vergence within the past few Gyr, although assigning a 
precise time to this phenomena is challenging due to un- 
certainties in AGB star modeling and the modest time 
resolution afforded by the data. The clearest differences 
between the morphological types can been seen in the 
most recent 1 Gyr, where the typical dSph, dl, dTrans 
and dSpiral formed - 2%, 8%, 4%, and 5% of their total 
stellar mass. 

The dwarf galaxies in the ANGST sample show a 
strong morphology-density relationship. This suggests 
that internal mechanisms, e.g., stellar feedback, cannot 
solely account for gas-loss in dwarf galaxies, as the corre- 
sponding models are unable to produce this observed re- 
lationship. Instead, we find qualit ative consistency wi th 
the m odel of 'tidal stirring' (e.s... iMaver et al.ll2001al lbl. 
I2006D . which can broadly explain the extended SFHs as 
well as the observed morphology-density relationship. 

A comparison with the cosmic SFH reveals gen- 
eral good agreement between the independently derived 
SFHs. The slower rate of SF for dwarf galaxies at in- 
termediate epochs may be tentative evidence of galaxy 
downsizing. However, broad uncertainties in extinction 
corrections and AGB star models may be able to explain 
the offset. Further, the mean measured SFHs are in- 
consistent with single valued exponential models of SF 
(i.e., r models), and may require more complex or multi- 
component models. 

We also identify 12 dTrans in the sample, based on the 
literature definition of present day gas fraction and SF 
as measured by Ha (e.g.. lMateoi ll998). Within this sam- 
ple of dTrans, we find that galaxies with high gas frac- 
tions are associated with more isolated galaxies, while 
those with lower gas fractions are less isolated. This 
suggests that there are two mechanisms that can pro- 
duce the observed dTrans characteristics: the isolated 
gas-rich galaxies are simply between episodes of SF due 
to the stochastic nature of SF in low mass galaxies, while 
the less isolated galaxies could be in the process of inter- 
acting with a more massive companion. 
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APPENDIX 

A. OPTIMIZING TIME RESOLUTION OF CMD-BASED SFHS 

Meaningful comparison of SFHs across the ANGST sample requires a uniform time binning scheme. This necessitates 
balancing the wealth of information about the ancient and intermediate epochs from the deepest CMDs (e.g., Antlia; 
50% completeness of Mj?8i4W ~ +1.75) with the course leverage of the shallower CMDs (e.g., UGC 4483; 50% 
completeness of Mpsuw ^ — 1-7; see Table [2|). As a compromise between the two extremes, we fixed a sample-wide 
time binning scheme using a galaxy with a 50% completeness value near the median of the sample (DDO 6; Mpgi4\Y 
~ 0). To test for the optimal time resolution, we constructed CMDs of single age stellar populations, convolved the 
simulated photometry with the artificial stars and photometric limits of DDO 6, and ran the SFH recovery program 
on each of these simulated CMDs. This process was conducted on single age CMDs of 0.5, 1.5, 2.5, 4.5, 8, and 12 Gyr 
in age. 

The purpose of this exercise is to provide a simple test for the robustness of the final time bins. While coarser time 
bins encompass a larger fraction of the input SF, there are generally fewer of them, which provides less information 
about pattens of SF over time. In this test, we deemed a good recovery as one in which the input SFH was within the 
uncertainties of the recovered SFH. In Figure I13[ we show the cumulative SFHs, where the input SF is indicated by 
the dashed magenta line and the recovered SFH by the solid black line. Uncertainties were computed using 50 Monte 
Carlo realizations as described in S}3l including identical shifts in Mboi and Te/ / • 

Initial tests suggested that six broad time bins of 0-1, 1-2, 2-3, 3-6, 6-10, 10-14 Gyr would provide suitable time 
resolution across the ANGST sample. As shown in Figure [131 the general agreement between the input and recovered 
SFHs is good and falls within the uncertainties. Clearly, the best recovery is seen in the 0-1 Gyr bin, where the young 
MS and massive core and blue helium burning stars provide secure leverage on recent SF. The oldest bin also shows 
similarly good results. 

The results in the intermediate bins are generally reasonable, although not as r eliable as the youngest and oldest 
bins. Much of this can be attributed to uncertainties in AGB star models (e.g., IGirardi et al.ri2010t iMarigo et al.l 
[20Tnh . Although, we used the same models to generate and recover the single age populations, photometric depth, 
systematic stellar model uncertainties, and blurring of distinct features due to observational effects all contribute to 
uncertainties in the recovered SFH (see Ap pendix [B I . AGB populations of different ages often have similar colors 



and magnitudes on optical CMDs (e.g.. [Gallart et al.ll2005HWeisz et al.ll2008f ). thus it can be difficult to confidently 



determine the epoch of SF at intermediate ages (1-10 Gyr ago). This limitation is clearly illustrated in Figure fT3l 
which shows broad uncertainties on the recovered SFHs during intermediate times. However, we still find that the 
input SFHs are generally within the recovered SFH error bars. Additionally, > 50% of the SF is recovered in the 
correct time interval. Although it would be possible to combine bins at intermediate times to increase the accuracy 
of the recovered SFRs, we believe that this scheme will permit future comparisons with SFHs derived using improved 
AGB star models (e.g.. IGirardi et al.|[2010l ). 

For this time binning scheme, we find that reported uncertainties in the SFRs typically decrease for CMDs deeper 
than DDO 6, and increase for shallower CMDs, in agreement with expectations. We thus conclude that this scheme 
is adequate for comparison of SFHs across the ANGST sample. 

B. EXPLORING UNCERTAINTIES IN CMD-BASED SFHS 

Photometric depth is an important consideration when interpreting the accuracy of a measured SFH. Varying 
photometric depths impact both the number of stars on a CMD and determine the presence of age sensitive CMD 
features (e.g., ancient MS turn-off, horizontal branch). Intuitively, a CMD with a brighter photometric limit has less 
information available than a deeper CMD, and the derived SFH is thus more uncertain. However, quantifying the 
precise impact of these uncertainties, such as the amplitude and ages affected, is challenging as varying photometric 
depths can amplify effects of uncertainties in the stellar models used to measure SFHs. In this section, we explore the 
effects of photometric depth on the accuracy of SFH recovery. We demonstrate the effects in two regimes: one where 
the only variable is photometric depth, and the other where we vary both photometric depth and stellar evolution 
models. 

We first analyze the accuracy of recovered SFHs as a function of only photometric depth. That is, we construct 
synthetic CMDs at select photometric depths, then attempt to recover the input SFH using identical parameters 
(e.g., IMF, binary fraction, filter combination) and, in this case, the same stellar evolution models. More concretely, 
we constructed CMDs at six different photometric depths {My = +4, +2, +1, 0, —1, —2), using a constant SFH 
(log(t) =7.4-10.15, with a time resolution of ~ 0.1 dex), a fixed metallicity, and the BaSTI stellar evolution models. 
Each CMD was populated with ~ 10^ stars in order to minimize the contribution of random uncertainties due to the 
number of stars used to measure the SFR in each time bin; For SFHs measured from observed CMDs, the uncertainties 
for Poisson sampling of the CMD are already well characterized by our Monte Carlo tests. We then recovered the SFH 
of each CMD using the BaSTI stellar evolution models, to ensure that the only variable being tested is photometric 
depth. 

In Figure [HI we compare the input (black lines) and recovered (colored lines) cumulative SFHs at each photometric 
depth. Overall, the recovered cumulative SFHs are in excellent agreement with the input SFHs at all photometric 
depths. The maximum deviation between the input and recovered SFHs at any photometric depth is ~ 4%, which is 
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Figure 13. A comparison between simulated input single age cumulative SFHs and the recovered SFHs (solid black line), 
demonstrating the time binning scheme for the ANGST sample. Here the magenta dashed line indicates the shape of a perfectly 
recovered cumulative SFH in the correct time bin. Delta function SFIfs were created at 0.5, 1.5, 2.5, 4.5, 8, and 12 Gyr ago 
and convolved with the artificial star tests of DDO 6, a galaxy representative of the median photometric depth in the ANGST 
sample. These single age CMDs were then recovered in a manner identical to the SFH recovery method, including 50 Monte 
Carlo realizations for error analysis. The grey dot-dashed line shows indicates 50% of the total stellar mass formed. We find 
general consistency between the input and recovered SFHs, particularly in the 0-1 and 10-14 Gyr time bins. Intermediate ages 
show a lower fraction of stars recovered in the desired bin, however all have recovery fraction > 50%. This broad time binning 
scheme allows us to balance a variety of photometric depths across the sample with temporal information about patterns of SF. 
See Appendix \K\ for a more detailed discussion. 



consistent vifith the expected Poisson precision of 1/ ^/lN) , where N is the number of stars used to measure the SFR 
in a given time bin. This exercise demonstrates that if all the underlying stellar models are known exactly, then the 
accuracy of the recovered SFH only depends on the number of stars in the CMD, and not the p hotometric depth. The 
same results are found when using different stellar evolution models, e.g., Padova, Dartmouth ([Potter et al.ll2008f) . to 
conduct this exercise. 

Systematic uncertainties are introduced into SFH measurements by uncertainties in the selected stellar models. The 
luminosity, color, and number density of specific CMD features provide leverage on the SFR at different epochs. How- 
ever, stellar models are not always self-consistent when trying to model multiple observed features (e .g., reproducing 
colors consistent with observations for both the horizontal branch and RGB; e.g.. iGallart et aLll2005[ ). The effect on 
a measured SFH is that the SFR may be systematically shifted into a particular time bin, depending on the stellar 
model used and photometric depth of the CMD (i.e., which particular age sensitive CMD features are available). We 
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refer the reader to iGallart et alj (|2005( ) for a comprehensive review on the adequacy of stellar evolution models for 
reproducing CMD features. 

To test for systematic effects on measured SFHs, one ideally wants to test for differences between the stellar models 
and 'truth', that is, stellar population characteristics based on the exact physics governing stellar evolution in nature. 
However, current stellar evolution models represent the best physical descriptions we have of nature. By measuring 
how different stellar models impact a measured SFH, we can get a sense of the amplitude of systematic uncertainties 
for SFH recovery. 

To examine the magnitude of systematic effects, we again construct six CMDs at selected photometric depths each 
containing ^ 10^ stars, assuming a constant SFH, fixed metallicity, and the BasTI stellar evolution models. However, 
the SFHs are now recovered with the Padova stellar models. Thus, in this case the differences between the input and 
recovered SFHs are indicative of the systematic effects introduced by choice in stellar model. 

In Figure [Tsl we see the differences in the input (black lines) and recovered (colored lines) mean cumulative SFHs are 
more substantial than when identical models were used. For the deepest CMD considered {My =+4), the agreement 
between the input and recovered CMD is excellent. Here, the ancient MS turnoff provides a reliable constraint on the 
ancient SFH, meaning the systematics in progressively younger time bins are also known more precisely. For shallower 
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Figure 14. A test of the effects of only photometric depth on the accuracy of SFH recovery. The black solid line is the input 
constant SFH, while the colored lines are the recovered SFHs at six select photometric depths. Synthetic CMDs were generated 
using the BaSTI stellar models, a constant SFH, and a fixed metallicity. The SFHs of the synthetic CMDs were then recovered 
using identical parameters, including the BaSTI stellar models. The input and recovered SFHs are in excellent agreement at 
all photometric depths. This implies that if the underlying stellar model is known exactly, then a SFH can be recovered to a 
precision limited by the number of stars in the CMD. 
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CMDs that do not include the ancient MS turn-off, we see significant deviations between the recovered and input 
SF Hs. These difference arise primarily in the treatment of such features of the horizontal branch, red clump, and RGB 
(see IGallart et al][2005l for a more detailed discussion). 

We quantify the magnitude of the discrepancies between the Padova and BaSTl models by examining the absolute 
value of the differences in input and recovered cumulative SFHs fPigure [TB)) . Of interest are the general order of 
magnitude variations, and not a specific point by point analysis. In the case of the shallowest CMDs (My =—2, — 1, 
0), the ancient SF is preferentially recovered at intermediate ages, with a typical magnitude in the difference of input 
and recovered SFH of ~ 20%. For the deeper CMDs (My =+1, +2), SF at the oldest times tends to be over-estimated 
by up to ~ 40%. At all photometric depths, the systematic effects within the last ~ 2 Gyr are < 10%. 

The exercises we have conducted demonstrate the importance of systematic uncertainties in SFH recovery. We 
caution that our results are actually measuring the systematic differences between the Padova and BaSTI models, and 
only serve a a proxy for the difference between a given model and 'truth', i.e., observed CMDs. That being said, the 
methodology of this type of analysis provides a framework for exploring the effects of systematic uncertainties. The 
primary limitation to this type of analysis is the number of models available that span the entire age/metallicity range 
needed to derive SFHs in nearby galaxies. As the number of models that sample a broader parameter space increases, 
it will be possible to gain more leverage on the effects of systematic uncertainties. 




Figure 15. A test of the effects of photometric depth and stellar model differences on the accuracy of SFH recovery. The blacfc 
solid line is the input constant SFH, while the colored lines are the recovered SFHs at six select photometric depths. Synthetic 
CMDs were generated using the BaSTI stellar models, a constant SFH, and a fixed metallicity. The SFHs of the synthetic 
CMDs were then recovered using identical parameters, and the Padova stellar evolution models. The input and recovered SFHs 
show excellent agreement for the deepest CMD. Discrepancies in the SFHs at shallower depths are indicative of the systematic 
effects introduced by the choice of stellar model. 
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Figure 16. An estimate of the amplitude of systematic uncertainties on SFH between the BaSTI and Padova stellar evolution models. Specifically, 
this figure highlights the absolute difference in the simulated input constant SFH (constructed with BaSTI) and the recovered SFH (solved with 
Padova) at six select photometric depths. In general, deep CMD including the ancient MS turn-off (A/v — +4) show very little systematic difference. 
SFHs from slightly shallower CMDs show SF at ancient times is generally overestimated by the Padova models, while the SFHs recovered from the 
shallowest CMDs place too much SF at intermediate epochs. We caution that these results demonstrate the differences only between these two 
stellar models, and may not quantitatively reflect the differences between a given model and 'truth', i.e., observed CMDs. 



